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ABSTRACT 



Aims. We analyze new multicolor light curves and recently published radial velocity curves for close binaries QX And, RW Com, 
MR Del, and BD +07" 3142 to determine the physical parameters of the components. 

Methods. The light curves are analyzed using a binary star model based on Roche geometry to fit the photometric observations. 
Spectroscopic parameters, such as the mass ratios and spectral types, were taken from recent spectroscopic studies of the systems in 
question. 

Results. Our findings provide consistent and reliable sets of stellar parameters for the four studied binary systems. Of particular 
interest is the BD +07° 3142 system, since this is the first analysis of its light curves. We find that it is an overcontact binary of W 
UMa type and W subtype, and that each component has a large cool spot in the polar region. QX And is an A subtype, and RW Com 
a W subtype W UMa binary, and in both systems we find a bright spot in the neck region between the components. MR Del is a 
detached binary with a complex light curve that we could model with two cool spots on the hotter component. 
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In this study we present simultaneous solutions of new, high- 
quality CCD light curves* and recently published radial velocity 
curves of four close binaries, QX And, RW Com, MR Del, and 
BD +07° 3142. While one of the systems, MR Del, is a close, 
detached system with a third light contribution from its visual 
companion, others have parameters consistent with the charac- 
teristics of the W UMa type systems. QX And belongs to the A 
subtype and RW Com to the W subtype. The analysis for BD 
+07° 3142, published here for the first time, shows that it is a W 
subtype W UMa system. 

The information about observations and a general descrip- 
tion of the light-curve analysis are given first, followed by the 
detailed history of previous study, the results of this work, and 
discussion for each star. 



2. Observations 

CCD photometric observations of the close binaries QX And, 
RW Com, MR Del, and BD +07" 3142 were obtained with 
the Apogee ALTA U47 CCD camera attached to the 40 cm 
Schmidt-Cassegrain telescope located at the Ankara University 
Observatory (AUG) by using BVR filters in close accordance 
with the Johnson-Cousins standard system. After calibrating all 



the images (bias, dark, and flat correction) and using the IRAF 
task CCDPROC, we performed aperture photometry in the usual 
manner with the relevant tasks in the aperture photometry pack- 
age (apphot) in IRAJjQ- Both BVR magnitudes and their errors 
were computed as variable minus comparison. The nightly ex- 
tinction coefficients for BVR magnitudes were determined from 
the observations of comparison stars. Orbital phases of the light 
curves were computed based on the light elements calculated us- 
ing the times of minima recently published in the literature. 

The log of the observations, information related to the com- 
parison stars and the check stars used in the observations, and 
the light elements used for the calculation of photometric phases 
of the light curves are given in Tables 1, 2 and 3, respectively. All 
the the times of minima co vered by our observations wer e calcu- 
lated using the method of iKwee & van Woerden I dl956l) . based 
on unweighted BVR observations. Times of minima being pub- 
lished for the first time are listed i n Table 4. For RW C om, we 
used times of minima calculated bv l Yilmaz et al. I (l2009t) . 



3. Light curve analysis 

We analyzed al l the light curves using the program by 
Diurasevi c et al. I dl992al) . generalized for the c ase of an 
overcontact configuration fDiurasevi c et al. I IT998). The pro- 



Send offprint requests to: G. Djurasevic 



'The data is available in electronic form at the CDS, via anonymous ftp 
at cdsarc.u-strasbg.fr or via http://cdsweb.u-strasbg.fr/Abstract.html 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 



for Research in Astronomy, Inc., 
National Science Foundation. 



under cooperative agreement with the 



1 



G. Djurasevic, M. Yilmaz, O. Ba§turk et al.: QX And, RW Com, MR Del and BD +07° 3142 



Table 1. The observation log. 



Table 4. New times of minima derived from our B VR observations. 



System 


Obs. date 


Nightly 


mean errors [mag^ 


and integration times 


[sec] 






B 




v 








QX And 


2009/08/17 


±0.005 


15 


±0.004 


5 


±0.004 


3 




2009/09/25 


±0.005 


15 


±0.004 


5 


±0.004 


3 




2009/09/29 


±0.005 


15 


±0.004 


5 


±0.004 


3 


rv y v V— -urn 


7009/09/0? 


±0.010 


40 


±0.005 


17 


±0.004 


10 




2009/03/70 


±0.010 


40 


±0.005 


18 


±0.004 


10 




2009/03/23 


±0.010 


35 


±0.005 


20 


±0.004 


10 




2009/04/08 


±0.010 


40 


±0.005 


17 


±0.004 


10 


MR Del 


2009/06/16 


±0.003 


10 


±0.003 


3 


±0.003 


2 




2009/07/29 


±0.003 


25 


±0.003 


4 


±0.003 


2 




2009/08/15 


±0.003 


25 


±0.003 


4 


±0.003 


2 




2009/09/25 


±0.003 


20 


±0.003 


8 


±0.003 


3 


BD 


2009/08/07 


±0.003 


30 


±0.002 


9 


±0.002 


5 


+07° 3142 


2009/08/09 


±0.003 


30 


±0.002 


9 


±0.002 


5 




2009/08/19 


±0.003 


30 


±0.002 


10 


±0.002 


5 




2009/08/31 


±0.003 


35 


±0.002 


10 


±0.002 


6 




2009/09/01 


±0.003 


35 


±0.002 


10 


±0.002 


6 



Variable 


Times of min. 


Type 


QX And 


2455061.4672(1) 


I 




2455101.4475(2) 






2455104.3342(2) 






2455104.5398(1) 


U 


MR Del 


2455042.5056(1) 


I 




2454999.4669(1) 


II 


BD +07° 3142 


2454688.3643(2) 


I 




2454686.2978(2) 


II 



Table 2. The light elements used in this study for each system. 





Epoch 








System 


(HJD+2400000) 


Period (days) 


Ref. 


QX And 


55101.4475(19) 


0.4121716(9) 


This 


study 


RWCom 


54272.8826(6) 


0.2373463 (1) 


This 


study 


MR Del 


55042.5055(7) 


0.5216903(1) 


This 


study 


BD +07" 3142 


54688.3631(12) 


0.275205(1) 


This 


study 



The fitting was done for a limited set of parameters, while 
the following parameters were fixed based on information from 
independent analysis of the radial velocity curves and plausible 
assumptions about the systems. 

- The mass ratio q = m c /mi, (where the subscripts h and 
c stand for the hotter and cooler component) for each 
system was fixed bas ed on the radial ve l ocity studies by 
iPribulla et all d2009bl) and lRucinski et al. I (120081) . 

- The temperature of the hotter component Th for each system 
was fixed according to its spectral typ e and the revised theo - 
retical T e ff-spectral type calibration by Marti ns et al. I (120051) . 

- The gravity-darkening coefficients /?/, ( . and the albedos Au, c 
for each system were fixed to their theoretical values accord- 
ing to the the temperature of the hotter component and a pre- 
liminary estimate of the temperature of the cooler compo- 

Table 3. The comparison and check stars used during the observation dPelfch close binary. 

- The ratio of the rotation rate to the Keplerian orbital rate 

fh,c = Wh,c/wn for all systems was fixed in accordance with 
the assumption that the rotation of the components is syn- 
chronous with the orbital revolution since the tidal effect 
are expected to contribute to the synchronization of the rota- 
tional and orbital periods. 

The limb darkening in the latest versio n of the mode l fol- 
lows the nonlinear approximation given bv lClaretl d2000h and 
the limb-darkening coefficients for the corresponding passbands 
were interpolated from Claret's tables according to the current 
values of T e ff and log(g) in each iteration. 

The results of the analysis are given for each star in a table 
(see Tables [51 16] [71 and [8]), where n is the total number of the 
B, V, and R observations; E(0 - C) 2 - the final sum of squares 
of residuals between observed (LCO) and synthetic (LCC) light- 
curves; <x - the standard deviation of the residuals; q = m c /mh 
- the mass ratio of the components; T niC - the temperatures of 
the hotter and the cooler component; /3h, c , Ah jC , /h,c - the gravity- 
darkening exponents, albedos and nonsynchronous rotation co- 
efficients of the components; A s , 6 S , A s , and <p s - the tempera- 
ture coefficients, angular dimensions, longitudes, and latitudes 
(in arc degrees) of the spots (if present in the model); Fh ;C - 
the filling factors for the critical Roche lobe of the hotter and 
cooler component; i [°] - the orbit inclination (in arc degrees); 
L3/(Lh +L C +L3) - the third light contribution; Qh,o> an d 
Q. out - the dimensionless surface potentials of the system compo- 
nents and of the inner and outer contact surfaces; / ver[%] - the 
degree of overcontact; Rh jC - the polar radii of the components in 
units of separation; Lh/(Lh + L c ) - the luminosity (B;V;R) of the 
hotter star; A4h,c[M ], ■Rh.JRo] - the stellar masses and mean 
radii of the components in solar units; log gh, c - the logarithm 
(base 10) of the system components effective gravity; M,' c , - the 





Variable 


Comparison 


Check 




QX And 


GSC 2816-2039 


GSC 2816-1820 


02000: 


01*57'"57'.78 


01*57 m 47 s .14 


0l''57"'31 s .94 


62000: 


+37°48'22".4 


+37°47'30".3 


+37°53'43".2 


Spec. Type: 


F5 


F0 


F3 


V(mag): 


11.43 


10.9 


11.2 




RW Com 


GSC 1991-1657 




ff 2000 : 


i2*33"w.28 


12*33«03».94 




62OOO: 


+26°42'58".4 


+26°37'27".6 




Spec. Typ.: 


G2 






V(mag): 


11.00 


11.86 






MR Del 


GSC 518-13 


GSC 518-447 


02000: 


20*31"'13 s .46 


20''31'"24 I .18 


20''3r23 s .53 


62000: 


+05°13'08".5 


+05°12'50".l 


+05° 18' 03". 1 


Spec. Type: 


K0 




F2 


V(mag): 


11.01 


9.86 


9.73 




BD +07" 3142 


GSC 380-167 


GSC 380-429 


02000: 


16''10'"03 s .21 


16' , 20" , 09 I .22 


16'' 20™ 12 s . 84 


62000: 


+07° 07' 28" .7 


+07°03' 54".6 


+07° ll' 17". 5 


Spec. Typ.: 


G5 


F8 




V(mag): 


9.89 


9.15 


10.6 



gram is ba sed on the Roche model a nd the principles de- 
scribed by IWilson & Devinnevl d!971l) . The system param- 
eters are estimated by applying an inverse-probl em method 
( Diurasevic et al. I fT992b). based on the modified Marquardt 
(U961 " algorithm. The underlying binary system model allows 
various system configurations from detached to overcontact, in- 
cluding active spot regions on the components. The latest version 
of the model and the so lving procedure is described in detail by 
iDjurasevic et al. I(l2008l) . 
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absolute bolometric magnitudes of the components; and a or b[R Q ] 
- the orbital semi-major axis in units of solar radius. 

Although the mass ratio for each system was fixed when 
solving the inverse problem to the value adopted from the rel- 
evant spectroscopic study, the influence of the mass ratio uncer- 
tainty was incorporated in estimating the final parameter errors. 
Namely, we made trial runs by fixing the mass ratio to its high- 
est and lowest possible values (chosen according to the reported 
uncertainty), and the deviations of parameter values obtained in 
this way were combined with the formal fitting errors (in the 
sense of adopting the highest error). Tables [5] [6] [7] and [8] list 
the parameter uncertainties estimated in this way. Still, the real 
parameter uncertainties may be greater than the tabulated errors 
since a significant contribution to the errors may come from the 
uncertainty of the effective temperature of the hotter star, which 
was in each case fixed on the basis of its spectral type. 

The distances to the four systems were computed based on 
the apparent magnitudes taken from the SIMBAD database and 
the computed absolute magnitudes, with corrections for the in- 
terstellar absorp tion. We used t he bolometric corrections from 
the tabulation o f Fl owerl dl996l) . according to the computed ef- 
fective temperatures of the components for each system. The in- 
terstellar extinction values (A v ) in the V passband were com- 
puted using the reddeni ng values estimated f rom the infrared 
dust emission maps of ISchlegel et all (Il998l) and by assum- 
ing an extinction to reddening ratio of A V /E(B - V) of 3.1. 
ISchlegel et al. I d 1998b refer to the total absorption. 

In the following sections we give the detailed history of pre- 
vious knowledge and the results of the present study for each 
system in turn. 

4. QX And 

QX And (GSC 02816-01950, NGC 752-235) is a short-period 
eclipsing binary, loca ted in the intermedia te-age open cluster 
NGC 7 52 (PC 1 363 inlRuprecht et al. I1988L age=2.0 + 0.2x 10 9 
vr from [ E)aniel et al. 1119941) . The system is also an x-ray source 
( Belloni & Verbunt 1 1 9961). The li ght v ariability of the system 
was first noted by iJohnsonl (Q953), and iTwaro g I (1 1 9 8 3l) pointed 
out that QX And was probably a binary based on its location in 
the color-magnitude diagram. 

The first photometri c and spectroscopic observations of QX 
And were obtained by Milone et al. ( 1995). They analyzed the 
light and radial v elocity curves simultaneou sly with the Wilson- 
Devinney code (Wilso n & Devinnevl 1197 lb . and obtained a 
model of an A-type W UMa system with a contact parameter 
of / = 0.21 + 0.11 that fit the observations best. The authors 
estimated the spectral type of the system as F3-5, determined 
the absolute parameters of the components (Mi = 1.1 8M and 
Mi = 0.24M o ), and found the distance to the binary and the clus- 
ter to be 381 + 17 pc. The radial velocity data of QX And was 
derived using cr oss-correlation funct ions. Unfortunately, the RV 
observations of Mi lone et al. I d 19951) were of rather poor qual- 
ity, resulting in only marginally useful orbital parameters. The 
authors obtained four light-curve solutions (the spotless, one- 
spot, two-spot, and three-spot cases) based on their finding that 
QX And displayed night-to-night variations in its light curves. 
Although the three-spot model fit the observations best, they de- 
rived the absolute parameters from the spotless-case solution. 
The authors assumed all spots were located on the hotter com- 
ponent 

iMilone et al. I(fl9 95) located a spot on one of the components 
although they performed their solution under the assumption 
of a radiative envelope. Noticing the mistake in this approach, 



Table 5. Results of the simultaneous (BVR) analysis of the QX 
And light curves obtained by solving the inverse problem for the 
Roche model with a bright spot on the more massive (hotter) 



component. 


Quantity 


n 


2689 


S(0 - C) 2 


0.1996 


cr 


0.0086 


q = m c /rai, 


0.306 ± 0.009 


n 


6440 




0.5 




0.08 


/« = h 


1.0 


A hs = T bs /T h 


1.02 ±0.01 




19.6 ± 1.0 


As 


357.3 ± 1.4 


f>bs 





To 


6420 ± 20 




1.031 ±0.001 


in 


54.6 ± 0.2 




2.41 ±0.02 




2.48 ±0.02, 2.29 ±0.02 


/over[%] 


35.2 ±2.0 


R b [D = 1] 


0.47 ± 0.01 


MD = 1] 


0.28 ±0.01 


L h /(L h + L c ) 


0.74[B]; 0.74[V]; 0.74[R] 


M h [M ] 


1.47 ±0.05 


M C [M ] 


0.45 ± 0.02 


%[R Q ] 


1.46 ±0.02 


K[Ro] 


0.88 ± 0.02 


loggh 


4.27 ± 0.03 


log gc 


4.20 ± 0.03 




3.49 ±0.03 




4.60 ± 0.03 


a orb[Ro] 


2.89 ± 0.04 


rf[pc] 


416 ±46 



Oian et al. ] (120071) find that QX And is a deep overcontact bi- 
nary system with a high degree of overcontact (f = 55.9%) and 
a low mass ratio (q = 0.2327). From the published times of min- 
ima, they point out that the orbital period is changing at a rate 
of dp/dt - 2.48 x 10~ 7 days/yr, which indicates a mass transfer 
from the less massive to the more massive component. Because 
the light curve they obtained displayed no asymmetry, they did 
not locate any spo t s on eit her of the components. 

IPribulla et al. (2009b) obtained radial velocity curves of 
both components and determined a mass ratio of q = 0.306. The 
authors estimated the spectral type of the system as F8. This re- 
ce nt mass ratio is more reliable than the one (q = 0.203) found 
bv lMilone et al. I d 1995b . The difference in the mass ratios found 
by these two authors is substantial, which is why new analysis, 
based on a reliable, spectroscopically estimated mass ratio and 
quality photometric observations was required for a better un- 
derstanding of this system. 

The results of our photometric analys is, based on the spectro - 
scopically estimated mass ratio given by Pribulla et al. I (12009b). 
are presented in Table [3] Figure Q] shows the observed (LCO) 
and the synthetic (LCC) light curves in the B, V, and R filters 
(upper left), the B - V and V - R color indices (lower left), the 
O - C residuals (upper right) and the geometrical model of the 
system in representative phases 0.3 and 0.8 (lower right). The 
parameters are given with the uncertainties estimated by com- 
bining the formal nonlinear, least-squared fitting errors with the 
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Fig. 1. Observed (LCO) and synthetic (LCC) light curves of QX And, obtained by simultaneous analysis of the B, V, and R obser- 
vations, A{B - V), A(V - R) color curve, final O-C residuals, and the view of the system model at the orbital phases 0.30 and 0.80, 
obtained with parameters from the inverse problem solution. 



errors arising from the uncertainty of the spectroscopic mass ra- 
tio (q = 0.306 + 0.009), as described in Section[3] 

Our solution shows that the less massive and slightly cooler 
component eclipses the more massive star in the primary min- 
imum (at phase zero), placing QX And in the A subclass of 
W UMa type systems, with a significant overcontact degree of 
/over ~ 35%, and a small temperature difference between the 
components. The distance estimation from the absolute mag- 
nitude computed with o ur model ag r ees wi th the estimation of 
380 + 36pc obtained by iBilir et al. I (120051) . using the parallax 
value. 

We did not observe any night-to-night variations in the light 
curves. The absence of significant phase dependence of color 
curves in Figure Q] indicates that there is good thermal contact 
between the components. However, for successful fitting of the 
slight asymmetry in the light curves, the optimal model of the 
system needs a bright spot located in the neck region of the more 
massive, hotter component. This result agrees with the analysis 



bv lOian et al. I d2007h . who find that mass transfer from the less 
massive to the more massive component, indicated by their O-C 
analysis, could justify a possible hot spot in the neck region of 
the more massive component. With the large difference in the 
masses of components, their equal temperatures suggest a sig- 
nificant energy exchange through the neck region. 

Because of the specific light curve shape with almost equal 
depths of the primary and secondary eclipses and a relatively 
small amplitude resulting from the low orbital inclination (z « 
55°), we also tested the hypothesis that the system is in a W 
subtype W UMa configuration. However, this model gave a sig- 
nificantly worst fit to the observations. 

In comparison with previous s tudies of QX And by 
Milo ne et al. I d 1995b and lOian etaTI d2007l) . who used the un- 
derestimated mas s ratio of q » 0.2 obtained from the ra- 
dial velocities of Milone et al. (1995;), our solution yields a 
significantly more massive secondary {AA Cf Milone = 

0.24 M , 

A1 C , Quian = 0.27 M B , and M Ci this study = 0.45 M ), an overcon- 
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tact parameter that is between their estimated values (fymoae — 
21%, /Q U j an = 55%, and / t hi s stu dy = 35%), somewhat larger radii 

Of the Components ("R^ Milone = 0.68 Rq, ??h, Milone — 1.39 R Q , 

and this study = 0.88 R Q , 7? h , this study = 1-46 R ), and a lower 
inclination (;' M iione ~ 58°, / Quian * 56°, and / this study « 55°). 
The temperatures, absolute bolometric magnitudes, and the esti- 
mated distance of the system agree well between these two stud- 
ies and our own. 



5. RW Com 

RW Com (HIP 61243) is a sho rt-period eclip sing binary. Its vari- 
ability was first recognized bv lJordaril d 1923b . who discovered it 
while observing the Ce pheid variable S Comae and classified it 
as B Lyrae type binary. IStruve I ([1950) performed the first spec- 
tral observations of the system and suggested that the emission 
components in the Ca II absorption line, observed to be strong 
in both conjunctions and weak in quadratures, originated in both 
star^ 

lO'ConneUI (1 19511) noted t he asymmetry and variability of 
the light curves, and iMilone I (Q976) interpreted them as ev- 
idence of mass loss from the system, primarily through the 
outer Lagrangian point, f acing the observer at primary mini- 
mum. Milone et al. Hl98Ch constructed light curves based on all 
data up to the time of their study, and also noted the variability 
in the light curves and the differences between the light levels of 
maxima. The asymmetry was re ported as increasing with wave- 
length by |DavidgeitaL| ([1983) • 

From their spectroscopic study, IMilone et al. I d 19851) found 
the mass ratio of 0.34 and classified the system as a W-subtype 
W UMa bi nary of spectr a l type G5-G8. Based on the results of 
this study, Mil one et al. I d 1987b analyzed UBV light curves of 
the system with the Wilson-Devinney code. Fixing the temper- 
ature to three different values (5400 K, 5600 K, and 5800 K) 
according to the estimated spectral type range (G5-G8), they 
tested different solutions, including both cold and hot photo- 
spheric spot configurations, and gave the absolute parameters of 
the solution that fit the observations best. 

ISrivastaval dl987T) published the first O-C curve of the sys- 
tem based on the previously published times of minima, and 
found a cyclic variation in the orbital period. Rely ing on the 
closel y spaced data obtained between 1967 and 1986, ISrivastaval 
(119871) suggested that a third body with an orbital period of 16 
years was causing the period to vary in time. lOianl d2002l) . how- 
ever, found that the variation caused by the third body had a 
13.3 year period, and was superimposed on a decreasing trend 
that the author explained b y angular momentum loss from the 
system. Yang & Liu I d2003l) considered both mass exchange and 
mass loss mechanisms that could cause the observed decrease in 
the period but found that the cyclic variation superimposed on 
the decreasing trend was quasi-periodic; hence, it could only be 
attributed to m agnetic ac t ivity. 

iPribulla&Rucinskil d2006l) and iRucinski et al. I d2007t) did 
not detect the signature of a third body in their studies dedi- 
cated to detecting faint third bodies in contact binary systems. 
In addition, broadening functions constructed for radial velocity 
analysis bv lPribulla et al . (2009a) did not show any evidence of 
a third companion in the system. They also found a much h igher 
mass ratio (0.471) than that found by Mil one et al. I d 1985b . on 
which all the light c urve analyse s publis hed up to that time were 
based. Moreover, iPribulla et al. I d2009al) determined a different 
spectra l type (K2 V) for the system, and conclude that the obser- 
vations lMilone et al71 d 19851) used were subject to errors because 



Table 6. Results of the simultaneous (BVR) analysis of the RW 
Com light curves obtained by solving the inverse problem for 
the Roche model with a bright spot on the less massive (hotter) 



component. 


Quantity 


n 


1309 


S(0 - C) 2 


0.3345 


cr 


0.0160 


q = mh/m c 


0.471 ±0.006 




4900 


A c =A h 


0.5 


Pc =A 


0.08 


fc = h 


1.0 


Abs = T bs /T h 


1.09 ±0.02 




14.2 ± 1.0 


As 


169.9 ±3.0 


fbs 





T c 


4720 ± 20 


Fc 


1.007 ±0.001 


in 


74.9 ±0.1 




2.80 ± 0.01 




2.82 ±0.01 ,2.54 ±0.01 


/over[%] 


6.1 ± 1.0 


RAD = 1] 


0.42 ± 0.01 


R h [D = 1] 


0.30 ±0.01 


L c /(L C +L h ) 


0.60[B]; 0.61[V]; 0.62[R] 


Mc[M Q ] 


0.80 ± 0.02 


Mh[M ] 


0.38 ± 0.02 


7?c[R ] 


0.77 ± 0.02 




0.54 ± 0.01 


log gc 


4.57 ± 0.03 


log gh 


4.54 ± 0.03 


M Li 


6.24 ± 0.05 




6.82 ± 0.05 


fl orb[Ro] 


1.70 ±0.06 


4pc] 


91 ±5 



of the heavy blends in lines used in their cross-correlation anal- 
ysis, and the long exposure times they employed in the observa- 
tions. 

The results of our photo metric analysis, based o n spectro- 
scopic elements estimated by IPribulla et al. I d200 9a) are given 
in Table[6] Figure[2]shows the observed (LCO) and the synthetic 
(LCC) light curves in the B, V, and R filters (upperleft), the B-V 
and V - B color indices (lower left), the O-C residuals (upper 
right), and the geometrical model of the system in representa- 
tive phases 0.20 and 0.75 (lower right). The tabulated parameter 
uncertainties were estimated by combining the formal nonlinear 
least-squared fitting errors with the errors arising from the un- 
certainty of the spectroscopic mass ratio (q = 0.471 + 0.006), as 
described in Section[3] 

Our solution shows that RW Com belongs to the W subclass 
of the W UMa type systems, i. e., the less massive, but hot- 
ter component, is eclipsed in the deeper primary minimum (at 
phase zero). The system is in a slightly overcontact configura- 
tion with / OV ei ~ 6%, and the temperature difference between 
the components is small (AT = T n - T c a 180 K). The optimal 
model of the system includes a bright spot in the neck region of 
the less massive, hotter component, enabling successful fitting 
of the slightly asymmetric light curves. This assumption is in 
accordance with mass exchange from the more massive primary 
to the secondary component that would cause the observed de- 
crease in orbital period. At the same time, the hot region in the 
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Fig. 2. Observed (LCO) and synthetic (LCC) light curves of RW Com, obtained by simultaneous analysis of the B, V, and R 
observations, A(B - V), A(V - R) color curve, final O-C residuals, and the view of the system model at the orbital phases 0.20 and 
0.75, obtained with parameters from the inverse problem solution. 



neck zone of the common envelope, located on the less massive 
secondary, can be interpreted as one consequence of an intensive 
energy transfer from the primary to the secondary. 

The distance of the system, estimated from the absolute mag- 
nitude obtained with our model, is within the uncertainty limits 
of the dist ance based on t he new Hipparcos parallax of 85 + 1 8pc, 
given bv Ivan Leeuw en (2007). We found it unnecessary to in- 
clude a third light in the sol ution, and that supports the find- 
ings of lPribullaetan (|2009a), who found no spectroscopic ev- 
idence of a third body in the syste m. Compariso n of ou r results 
with the parameters reported by Mil one et al. I (1 19871) . which 
were based on the u nderestimated value of the mass ratio from 
iMilone etal. 1 (119851) . show that our solution gives higher masses 
of the components (M c , Milone = 0.56 M , Mu, Milone = 0.20 M 
and M c , this study = 0.80 M , M h , this study = 0.38 M G ), a smaller 
overcontact factor {fMiione = 17% and /,/„■, stu d y = 6%) and 
larger radii M iione = 0.67 R Q , % M iione = 0.41 R and 
■^c, this study = 0.77 R Q , t his study = 0.54 R ), with approx- 



imately the same inclination and considerably lower tempera- 
tures, which is a conseq uence of the cooler sp ectral type as- 
signed to this system by Pri bulla et al. (2009a). However, the 
absolute bolometric magnitudes and the estimated distances are 
in relatively good agreement. 



6. MR Del 

MR Del (GSC 00518-1755, HD 195434, ADS 13940) is the 
brighter comp onent of a visual bin ary system. It is an eclips- 
ing binary dCutispoto et al. I Il997h. separated from its visual 
companion by 1 .8 at 71° dMason et al. 11200 lh . The whole sys- 
tem w as classified as a halo population object or an old disk 
star by Cutisp oto et al. ( 1997 ) , based on its proper motion an d 
low metallicity (Carney et al. 1994; Sandage & Kowal 1986). 
This makes it an interesting o bject, because it displays soft X- 



ray emission ( Pye et al. I [T 9951) and photospheric spot activity 
(Cutispoto et al. Ifl997h despite its age. 
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Fig. 3. Observed (LCO) and synthetic (LCC) light curves of MR Del, obtained by simultaneous analysis of the B, V, and R obser- 
vations, A(B - V), A(V - R) color curve, final O-C residuals, and the view of the system model at the orbital phases 0.30 and 0.70, 
obtained with parameters from the inverse problem solution. 



An out-of-eclipse variation with an amplitude of 0.04 mag, 
observed in the V band light curve of the system, was interpreted 
as comi ng from a large spot c overing 23% of the primary star's 
surface (Cutis poto et al. 1119971) . Although the system also dis- 
played variations in X-ray emission, which were thought to be 
coming from the eclipsing compon ent of the visual binary , it was 
not correlated with the light curve (Cutispot o et al. |[l997l) . 

The eclipsing component of the visua l binary was later 
named MR Del by Kazarov ets et al. I (1 19991) . The dispute over 
the magnitude s of each of the com ponents of the visual binary 
was ended by iMason et al. 1 (1200 lh. who gave 9 m .49 for the V 
band magnitude of MR Del. Claus en et al. 1 (12001) obtained ubvy 
light curves of the system, and observed a flare event lasting for 
25 minutes. They also point out the night-to-night differences in 
the light curves, increasing in strength from y to u band, which 
they attribute to surface activity. 

Pribulla et al. I d2009bl) observed the system spectroscopi- 
cally. They computed the light contribution of the visual com- 



panion as L3/(L; + L2) = 0.51. The broadening functions they 
employed to obtain radial velocities did not show any evidence 
of photospheric spots. They found the mass ratio of 0.915 and 
emphasized the need for a new photometric analysis of the sys- 
tem to attain the absolute physical parameters. 

IZasche et al. I J2009) updated the light elements after having 
analyzed all photometric and astrometric data available for the 
system. According to the observational indicators, MR Del has 
properties similar to stars of BY Dra type or of short-period RS 
CVn stars. 

The results of our photo metric analy s is, ba sed on updated 
spectroscopic elements of Pribulla et al. (2009b), are given in 
Table [7] Figure [3] shows the observed (LCO) and the synthetic 
(LCC) light curves in the B, V, and R filters (upper left), the 
B - V and V - B color indices (lower left), the O-C residuals 
(upper right) and the geometrical model of the system in repre- 
sentative phases 0.3 and 0.7 (lower right). Table|7]lists parameter 
uncertainties estimated by combining the formal nonlinear least- 
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Table 7. Results of the simultaneous (BVR) analysis of the MR 
Del light curves obtained by solving the inverse problem for the 
Roche model with two cool spots on the less massive (hotter) 
component. 

Quantity 



•C) 2 



n 

E(0 
tr 

q = m^/m c 
T h 

A c =A h 

ft =A, 
fc = A 

Asi =T S1 /T h 

e s i 

-^S 1 
jPSl 

As2 = Ts2/Th 

f?S2 
>fS2 

T 
F b 
Fc 
HI 

L 3 /(L h + L c + L 3 )(B) 
L 3 /(L h + L c + L 3 )(V) 
L 3 /(L h + L c + L 3 )(R) 

ci b 

n c 

^ h [o = i] 

R C [D = 1] 
L h /(L h +L C + L 3 ) 
M C [M ] 
Alt, [Me] 

^c[R Q ] 
%[R ] 

loggc 
log gh 

K* 

<! 
£I 01 b[Ro] 

d[pc] 



2574 

0.3120 

0.0110 

0.915 ±0.012 

4900 

0.5 

0.08 

1.0 

0.70 ± 0.03 
19.7+2.2 
272.1 ± 1.2 
61.5 + 5.7 
0.70 ± 0.03 
14.8 ±0.7 
100.1 ± 1.3 
-27.2 ±3.8 
4400 ± 20 
0.74 ± 0.01 
0.62 ± 0.01 
78.7 ±0.1 
0.23 ± 0.03 
0.25 ± 0.03 
0.25 ± 0.02 
4.58 ± 0.09 
5.29 ± 0.06 
0.27 ±0.01 
0.21 ±0.01 

0.59[B]; 0.56[V]; 0.56[R] 

0.69 ± 0.07 

0.63 ± 0.06 

0.83 ± 0.02 

0.65 ± 0.02 

4.44 ± 0.02 

4.61 ±0.02 

5.92 ± 0.08 

6.92 ± 0.08 

2.99 ± 0.03 

63 ±9 



squared fitting errors with the errors arising from the uncertainty 
of the spectroscopic mass ratio (q = 0.915 + 0.012), as described 
in Section [3] 

Our model includes two cool spots on the more-massive, hot- 
ter component. The spotted model is supported by the X-ray ob- 
ser vations. Anoth e r activ ity indicator is the flare event observed 
by IClausen et al. which was most pronounced in the u 

band. In addition, there are night-to-night differences in the light 
curves, increasing in strength from the y to the u band, so cool 
spots can be expected on one or both components; however, the 
uniqueness of the spot locations obtained in our solution is ques- 
tionable to some degree. A good fit could not be obtained with 
the spots located on the less massive secondary. A model with 
spots on the primary fits the observations very well, but fits of 
similar quality can be obtained with different spot locations, and 
we could not decide definitely which of the possible solutions is 
the most appropriate. 

Having in mind the vicinity of the third visual component, 
which cannot be resolved spatially, we included a third light 



in the model. The estimated third light contribution to the to- 
tal light of the system (L3/(Lh + L c +L3) as 0.23) is signif- 
icant, but somewhat l ower than the contribution reported by 
iPribulla et al.l (l2009bl) (L 3 /(L! + L 2 ) = 0.51, which equals 
L 3 /(Lh +L C +L3) as 0.33 for the luminosities of the eclipsing 
components in our model). There is also a considerable differ- 
ence between our estimation of the distance of MR Del (63 +9pc) 
and the one computed based on the new Hipparcos parallax 
(49 + 6pc). The greater distance implies that the total luminos- 
ity of the system is overestimated in our model, which may be 
a result of complications introduced by the third light. We made 
tri al runs with the L3 valu e fixed to the spectroscopic estimation 
of IPribulla et al. I (l2009bl) and obtained plausible solutions, but 

the fits to the observations were sig nificantly worse. 

In comparison to the results that Cutis poto et al. I ( 1 19971) ob- 
tained from the analysis of their UBV(RI) and soft X-ray light 
curves, our study gives somewhat lower masses of the compo- 
nents (Ale, cutispoto = 0.64 M Q , M h , cutispoto = 0.72 M Q and 
Ale, this study = 0.69 M , M h , this study = 0.63 M ), different radii 
(%:, Cutispoto = 0.70 R s , 7?h, Cutispoto = 0.77 R Q and He tM S stu d y = 
0.83 R Q , 7?h, this study = 0.65 R e ), and a slightly lower inclination 
0'cutispoto ~ 80° and i Us study a 79°). We note that their analysis 
was only based on photometric data. 

7. BD +07" 3142 

BD +07° 3142 (GSC 00380-00247, ASAS 16200 3+0707.4 
PPM 162491)) was found to be a variable star by iPoimanskil 
(120021) during their analyses of the ASAS data. Since then, 
only one time of prim ary minimum has been published by 
Parim ucha et al. I d2007l). 

Rucinski et al. I (l2008l) observed the star spectroscopically to 



determine its radial velocity variation. They found a third com- 
ponent with a light contribution of L3/(Li + L2) = 0.50. Because 
of the proximity of the radial velocity of the third component to 
that of the center of mass of the binary system, they argue that it 
appears to be physically attached to the eclipsing pair. From the 
analysis of broadening functions, they determined the mass ratio 
of 0.662 and spectral type as K2V. They classified the system as 
a W UMa type binary. 

The results of our photometric analy sis, based on spectro- 
scopic elements of Rucinski et al. I (120081) . are given in Table [8] 
This is the first ever published photometric solution of this sys- 
tem's light curves. Figure [4] shows the observed (LCO) and the 
synthetic (LCC) light curves in the B, V, and R filters (upper 
left), the B - V and V - B color indices (lower left), the O - C 
residuals (upper right) and the geometrical model of the system 
in representative phases 0.2 and 0.8 (lower right). As in the pre- 
vious sections, the uncertainties given in Table|5]were estimated 
by combining the formal nonlinear least-squared fitting errors 
with the errors arising from the uncertainty of the spectroscopic 
mass ratio (q = 0.662 + 0.008), as described in Section[3] 

According to our analysis, BD + 7°3142 is an overcontact 
binary (/ 0V er ~ 10%) in which the hotter, less massive compo- 
nent is eclipsed in the deeper primary minimum, so the system 
belongs to the W subclass of W UMa type binaries. The light 
curve asymmetry can be explained by the presence of two cool 
spots, located near the polar regions of each component. 

Since the orbital period of the system is short (P w O^S), 
its components are in fast rotation. Owing to the tidal effects, 
we expect a synchro nization of the star rot a tion w ith the orbital 
period. According to lSchiisler & Solanski I (119921) . the presence 
of spots at high latitudes in such a situation can be explained by 
the dynamo mechanism for rapid rotators. 
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Fig. 4. Observed (LCO) and synthetic (LCC) light curves of BD + 7°3142, obtained by simultaneous analysis of the B, V, and R 
observations, A(B - V), A(V - R) color curve, final O-C residuals, and the view of the system model at the orbital phases 0.20 and 
0.80, obtained with parameters from the inverse problem solution. 



The inclination of the orbit was estimated to be ; « 72°, 
which means that both the primary and the secondary minima 
are due to partial eclipses. The third component (probably phys- 
ically attached to the eclipsing pair) contributes approximately 
48% to the total light. Such a strong third light effect has a seri- 
ous impact on the estimation of the basic system parameters and 
is particularly important for estimating the degree of overcontact 
and the orbital inclination. 

8. Conclusions 

We estimated the physical parameters of four binary systems - 
QX And, RW Com, MR Del, and BD + 7°3142 - from the anal- 
ysis of new, high quality CCD light curves in the B, V, and R fil- 
ters, and from the results of recent high-resolution spectroscopic 
studies. This is the summary of our findings. 

1 . QX And is a close binary belonging to the A subclass of W 
UMa type systems, in a deep overcontact configuration with 



the degree of overcontact of / a 35% and almost equal tem- 
peratures of the components (7/, = 6440-7T and T c = 6A22K). 
The inclination of the orbit is estimated to be i = 54.6 de- 
grees, and the distance of the system is d = 416 pc. The 
best-fitting model contains a bright spot in the neck region 
of the hotter (more massive) component. 

2. RW Com is a close binary of W UMa type and W subtype in 
a shallow overcontact configuration with the degree of over- 
contact of / « 6% and almost equal temperatures of the 
components (T h = 4900/T and T c = 4719/Q. The inclina- 
tion of the orbit is estimated to be i = 74.9 degrees, and the 
distance of the system is d = 91 pc. The best-fitting model 
contains a bright spot in the neck region of the hotter (less 
massive) component. 

3. MR Del is a detached close binary and a component of a vi- 
sual pair. Its complicated light curve was modeled with two 
cool spots on the hotter component. The estimated tempera- 
tures of the components are T h = 4900/T and T c = 43967T, 
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Table 8. Results of the simultaneous (BVR) analysis of the 
BD + 7°3142 light curves obtained by solving the inverse prob- 
lem for the Roche model with two cool spots, one on each com- 
ponent. 

Quantity 



E(0 - C) 2 
cr 

q = m h jm c 
T h 

A h =A c 

fa =A 
A = A 

A s i = T sl /T c 

#sl 
(Psl 

A S 2 = T s2 /T h 

9s2 
A S 2 
<Ps2 
T c 

F c 
HI 

L 3 /(L h + L c + L 3 )(B) 
L 3 /(L h + L c + L 3 )(V) 
L 3 /(L h +L c +L 3 )(R) 
f2 h , c 

*^in) ^out 
,/over[%] 

RdD = 1] 
R h [D = 1] 
L h /(L h + L c + L 3 ) 
At_[M Q ] 
Mh[M G ] 

%[R Q ] 

log go 

loggh 

ML 

Ka 

a rb[Ro] 

<flpc] 



1719 

0.4572 
0.0163 

0.662 ± 0.008 

4900 

0.5 

0.08 

1.0 

0.70 ± 0.02 
33.3+0.4 
190.4 ± 0.9 

53.8 ±0.7 
0.81 ±0.02 
33.1 ±0.5 
326.7 ± 1.8 

56.9 ± 1.0 
4640 ± 30 
1.014 ± 0.002 
71.6 ±0.3 
0.451 ± 0.004 
0.478 ± 0.003 
0.480 ± 0.004 
3.140 ±0.015 

3.18 ±0.02, 2.79 ±0.01 
9.5 ± 1.3 
0.396 ± 0.002 
0.328 ± 0.002 

0.274[B]; 0.265[V]; 0.268[R] 

0.74 ± 0.05 

0.49 ± 0.04 

0.81 ±0.01 

0.67 ± 0.01 

4.49 ± 0.03 

4.48 ± 0.03 

6.20 ± 0.09 

6.37 ± 0.09 

1.91 ±0.02 

76 ±6 



and the inclination of the system is ; = 78.4 degrees. Our 
estimation of the distance of the system is d = 63 pc. We 
find that the contribution of the third light in V band is 
L 3 /(L h +L c +L 3 ) = 0.23. 
4. BD + 7°3142 is a close, interacting binary. This is the first 
published analysis of its light curves. We find that it is an 
overcontact binary with a degree of overcontact of / « 10%, 
and we classify it as a W type W UMa system. The esti- 
mated temperatures of the components are T/, = 49007T and 
T c = A6A5K and the inclination of the system is ;' = 71.8 
degrees, with a considerable contribution of the third light in 
the V band of L 3 /(L h + L c + L 3 ) = 0.48. We estimate that the 
distance of the system is 76 pc. The best-fitting model con- 
tains two large cool spots, one on each of the components. 

Since all four binaries appear to display surface activity, long- 
term observational campaigns along with the application of 
imaging techniques to affirm the existence and fix the locations 
of stellar spots would make valuable contributions to our knowl- 
edge about these systems. The question of whether BD + 7°3 142 



belongs to a wide, visual system remains open and calls for fur- 
ther study. 
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